Laser-induced effects on the electronic features of graphene nanoribbons 
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We study the interplay between lateral confinement and photon-induced processes on the elec- 
tronic properties of illuminated graphene nanoribbons. We find that by tuning the device setup 
(edges geometries, ribbon width and polarization direction), a laser with frequency Q may either 
not affect the electronic structure, or induce bandgaps or depletions at ±/i,f2/2, and/or at other 
energies not commensurate with half the photon energy. Similar features are also observed in the 
dc conductance, suggesting the use of the polarization direction to switch on and off the graphene 
device. Our results could guide the design of novel types of optoelectronic nano-devices. 
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The extraordinary properties of graphene [IHSl led to 
an unprecedented narrovifing in the expected gap between 
the understanding of new phenomena and the develop- 
ment of disruptive applications [ij. Though originally 
focused mainly on pure electronic, mechanical or optical 
properties, much attention is now devoted to the inter- 
play between these variables [3] . Graphene optoelectron- 
ics P-[l0|. in particular, is one of the most active and 
promising fields with flagship applications including en- 
ergy harvesting devices [6j and novel plasmonic proper- 
ties!, El. 

Recently, the captivating possibility of controlling the 
electronic properties of graphene through simple illumi- 
nation with a laser field jlj, has been re-examined 
through atomistic calculations [ij], calculations of the 
optical response [H, [l^ and proposals for tuning the 
topolo gical pr operties of the underlying photon-induced 
states |l7l - l20j , among other interesting issues (2l| - |26j . 
The basic idea is that laser illumination may couple states 
on each side of the charge neutrality point inducing a 
bandgap at ±hfl/2, if the field intensity and frequency 
are appropriately tuned. This non-adiabatic and non- 
perturbative effect relies crucially on the low dimension- 
ality and peculiar electronic structure of graphene and 
has attracted much recent attention |17H20l |. Notwith- 
standing, most of these predictions were restricted to 
bulk graphene. One may wonder about the possible con- 
sequences of reduced dimensionality and quantum con- 
finement. 

Here, we address the effects of laser illumination on 
graphene nanoribbons and show that lateral confinement 
plays a crucial role: tuning the sample size and the 
direction between the laser polarization relative to the 
sample edges, linearly polarized light may or not induce 
bandgaps or depletions in the density of states and the 
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conductance spectra. Strikingly, for finite size samples 
these features may appear at energies different from in- 
teger multiples of ±fir2/2. This is in stark contrast with 
bulk graphene where the electrical response is insensitive 
to the linear polarization direction. Our results fill the 
gap in the understanding of the laser-induced control of 
the electrical response and may guide the design of new 
experiments on optoelectronic devices. 

Hamiltonian model and solution scheme. We consider 
an infinite graphene nanoribbon illuminated by a laser 
only in a finite region of length L and perpendicular to 
it (as shown schematically in Fig. [T]-a) . Electrons in the 
graphene ribbon are modeled through a nearest neigh- 
bours TT-orbitals Hamiltonian [1, H^: He — Yl,i^i^t'^i ~ 
j) [li,j'^t'^j +H.C.], where and Cj are the electronic 
creation and annihilation operators at site i, Ei are the 
on-site energies and the nearest-neighbors carbon- 
carbon hoppings amplitudes which are taken equal to 
7o = 2.7 eV Radiation is described through a time- 
dependent electric field E. By choosing a gauge such 
that E = —dA/dt, where A is the vector potential, the 
hopping matrix elements acquire a time-dependent phase 

according to: ^ij — 70 exp (^i^ J^' A{t) ■ dr^ , where $0 
is the magnetic flux quantum. 

Retaining non-perturbative and non-adiabatic correc- 
tions to the electrical response is crucial for the results 
presented hereafter. In this regime, Floquet theory pro- 
vides an appropriate framework. An efficient solution us- 
ing this scheme is used to obtain the average density of 
states and the dc component of the conductance, which 
is computed from the inelastic transmission probabili- 
ties in Floquet space [2^. The interested reader may 
find further generalities of the method in Refs. [28l . l29j |. 
while more technical details will be published elsewhere 
(30| . For a periodic modulation of the hoppings, the 
spectral and transport properties can be derived from 
the so-called Floquet Hamiltonian: Hp = — iHd/dt. 
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FIG. 1. (Color online) a- Scheme of the considered setup: 
A graphene nanoribbon illuminated by a mid-infrared laser 
within a finite region of length L. The unit cell for the case 
of an armchair ribbon is marked with a rectangle (dashed 
line) and contains 4A'^ atoms in the notation used in the 
text, h- Floquet quasi-energy dispersion for a small ribbon 
iy'N = 2). For the sake of better visualization, we show the 
quasi-energy spectrum without radiation (this allows to bet- 
ter distinguish the crossings between levels) and choose ar- 
bitrary frequency and intensity values. More realistic results 
are shown in Figs. [5] and O The scheme on the right shows 
the electronic states for two of such transitions. 



Such Hamiltonian has a time- independent representation 
in the Floquet space, which is the direct product between 
the usual Hilbert space and the space of time-periodic 
functions with the same period as the Hamiltonian Hp 
[28l Isij . Therefore, on top of the k label, our states 
have a second label which indicates the number n of 
photon excitations in the system. In the absence of ra- 
diation, the quasi-energies spectrum of Hp are given by 
e(k, n) — eo(k) + nMl (£o(k) is the spectrum of H^). 

Electronic properties of irradiated graphene nanorib- 
bons. The first question we address is whereas the re- 
sponse of graphene nanoribbons to a laser is sensitive 
to the (linear) polarization direction. While in the bulk 
limit both the conductance and the density of states are 
independent on the polarization direction, the picture 
turns out to change radically in confined geometries. 

Figure ^ shows the average density of states as a func- 
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FIG. 2. (Color online) Average density of states for armchair 
graphene nanoribbons with A'^ = 130 {a,b) and N — 129 (c,d) 
for two values of the laser power P: lmW//im^ (solid black 
line with orange shaded area) and lOmW/zim'^ (dotted red 
line). Panels a and c {b and d) are for polarization along the 
X direction {y direction). The energy corresponding to 
is marked with an arrow for reference in the top of the panel. 
The DOS in absence of radiation is also shown for comparison 
(dashed lines, grey shaded area). 



tion of the Fermi energy for a frequency in the mid- 
infrared regime {hfl = 140meV). Two different armchair 
nanoribbons sizes and polarizations are chosen: N = 130 
(a and 6), N = 129 (c and d); and linear polarization 
along the x (a and c) and y {b and d) directions. For 
N = 130 one sees the appearance of strong depletions at 
zLhU/2. These depletions are located at the same energy 
as the ones for the bulk system ^T3] and correspond to 
the excitation of an electron between the conjugate states 
at zLhil/2. In striking contrast, the DOS is restored at 
±hfl/2 for y-polarized laser, whereas new features occur 
at energies incommensurate with hfl. 

The DOS for the ribbon with N ^ 129 (c and d) also 
exhibit a laser-induced fragmentation of the spectrum, 
with in this case the observation of fully depleted en- 
ergy regions (bandgap) at ±hQ/2. One observes that by 
increasing the laser intensity, some DOS modifications 
are further enhanced (see for instance the DOS deple- 
tion around hil/2 for Fig. ^a), complemented by the 
emergence of new fine structure. Figure [3] shows the dc 
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FIG. 3. (Color online) DC conductance for armchair graphene 
nanoribbons with — 130 (a, 6) and A'' — 129 {c,d). Panels a 
and c (6 and d) are for polarization along the x direction {y 
direction) . 



conductance as a function of the Fermi energy position 
for the same cases shown in Fig. [21 Here, we see that 
the depletions in the DOS have the same conductance 
fingerprints. One can see that switching the polarization 
direction may produce a marked on- off ratio if the Fermi 
energy is appropriately tuned. 

To rationalize these differences, it is instructive to 
write He in a basis of independent transversal channels 
or modes as discussed in Refs. [H, [s^l- In the absence of 
radiation, the Floquet spectrum of the system is just the 
sum of the contributions from each of the modes £o(k) 
plus their Floquet replicas: e(k, n) = eo(k) + nhil. An 
illustration for a very small system is shown in Fig. [TJ-6, 
each mode contains an electron and a hole branch. At the 
crossings between the different lines, one may expect for 
the effects of radiation to be stronger (if it provides the 
necessary coupling between the corresponding states). 

From the Floquet spectrum in Fig. [TJ-6 one can see that 
there are two kind of crossings: the ones that connect an 
electronic (hole) state with a hole (electron) state belong- 
ing to the same mode plus or minus an integer number of 
photons (like the one marked with an open circle); and 
the ones that connect states in different modes (as the 
one marked with an open square in Fig. [TJ-&). In such 



cases, a non-vanishing matrix element of the Floquet 
Hamiltonian will give intra and inter-mode transitions 
respectively. Given the electron-hole symmetry of the 
spectrum, intra-mode transitions always connect states 
which are symmetric relative to the Dirac point, i.e. in- 
teger multiples of ±fi,$l/2. On the other hand, inter-mode 
transitions always couple states which are not symmet- 
rically located from the Dirac point, as can be seen on 
Fig. [TJ-6. A scheme showing these two type of transitions 
is shown in Fig. [TJ-6. 

For armchair graphene nanoribbons, it turns out that 
a laser with linear polarization along the transport direc- 
tion (x) does not mix these transversal modes, leading to 
features in the density of states only at ±nhn/2 as can 
be seen in Fig. [TJ-&. On the other hand, calculation of 
the matrix elements shows that when the polarization is 
along y, inter-mode processes are allowed and intra-mode 
ones are forbidden. Depletions or gaps appear now at the 
crossing between Floquet states corresponding to differ- 
ent modes leading to the features observed in Fig. 
and c. 

A complementary approach to this problem is possi- 
ble by using the k.p model, which could be accurate 
enough for medium-sized ribbons [33 |. A careful anal- 
ysis shows consistent results: inter-mode processes lead 
to gaps/depletions located away from ±hn/2 while po- 
larization along y suppresses the depletions at ±hfl/2. 
In the bulk limit, as the energy difference between sub- 
bands gets smaller, the crossings between Floquet states 
accumulate close to ±hft/2 leading to the same behav- 
ior for both polarizations (along x and y). A flavor of 
this can be seen in the red dotted lines in Fig. and 
c: The two depletions seen in Fig. [^J-d black line merge 
when increasing the laser power. 

Another interesting feature is that the metallic 
modes/subbands in armchair ribbons are quite insensi- 
tive to the radiation (as seen in Figs. [2]and|31). Hence, 
very small metallic armchair ribbons containing only one 
transport channel within the energy range of interest will 
not experience relevant changes in their electronic prop- 
erties. We emphasize that this is a peculiar property of 
armchair graphene nanoribbons (soj . 

Figures [5] and [3] correspond to a laser frequency in the 
mid-infrared, which gives an optimum playground to test 
these predictions in the laboratory. Going to higher fre- 
quencies may help to achieve device miniaturization but 
the size of the gaps and depletions at constant laser power 
diminishes, whereas at very low frequencies (THz) the 
gaps further decrease. Experiments would require tem- 
peratures below 20 K for P ~ 1 mW//im^. 

In summary, we show that the interplay between 
photon-induced inelastic processes and lateral confine- 
ment in graphene nanoribbons leads to diverse modifi- 
cations in the band structure and transport properties 
not evident in the bulk limit. In the case of moderate 
sized nanoribbons (ca. lOnni), the careful tuning of the 
polarization direction may widen the opportunities for 
achieving control of the electrical response in optoelec- 
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tronic devices. 
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